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We survey cooperative phenomena of a cavity system in which many atoms interact with a single 
quantized photon mode driven by an AC external field in a dissipative environment. It has been 
. . . known that an external field causes the optical bistability which is a non-equilibrium phase transition 

' for the balance of excitation and dissipation. On the other hand, a strong interaction causes the 

Dicke transition which takes place with a spontaneous appearance of photons and excited atoms in 
, the equilibrium system. We study the phenomena in full range of the strengths of the interaction and 

the external field. When the interaction is strong, the simple Lindblad equation cannot reproduce 
' the correct equilibrium state. In order to treat the highly interacting region, we provide an extended 

, master equation which realizes the ground state in the case of no external field. We study the Dicke 

■ model and the Tavis-Cummings model under an AC external field by using the master equation and 

obtain the two following results: first, we present phase diagrams of the stationary state of the two 
models and we find that the rotating wave approximation causes a qualitative different nature of 
the phase diagrams when the interaction and the external field are strong. Second, we find a novel 
symmetry breaking phenomenon induced by the AC external field in the Dicke model. 
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§ ■ I. INTRODUCTION 

^3 ■ 

^ '- The effect of the interaction between photons and atoms has been studied for a long time. A cavity confines 
photons in a finite region with mirrors, which enhances the interaction between photons and atoms. When a single 
' atom is regarded to be confined in a cavity, the effect of interaction has been studied by using the Jaynes-Cummings 
model [l|, , where a single quantized mode of the cavity photons couples with a two-level system. For the system with 
' many atoms, the extension of the Jaynes-Cummings model has been introduced as the Tavis-Cummings model 3|. 
With increasing the number of atoms, the hybridization of the two-level systems and photon system is enhanced 4|. 
In experiments, we observe the effect of interaction as cavity ringing phenomena [a and the vacuum field Rabi 
splitting 0, 0]. Recently, this coupling leads to the ultra-strong coupling regime [8, Q, and also attracts much 
attention as a possible method to control to exchange quantum information between two-level systems and photon 
^ ; system [iMl- 

We also observe cooperative phenomena due to the interaction in a cavity. Even if there is no direct interaction 
between atoms, an effective long-range interaction appears among atoms via a common mode of cavity photons. It 
has been pointed out that the system shows two kinds of phase transitions: the Dicke transition and the optical 
bistability. The former is regarded as a ground-state phase transition as a function of the coupling strength at zero 
^ , temperature When the coupling strength is weak, the ground state of the system is a state without excitations 

^ ' of photons and atoms. When the interaction exceeds a critical value, photons and dipole moments of atoms appear 
spontaneously as a continuous transition. The latter is a kind of non-equilibrium phase transition under AC external 
field in contact with thermal baths. It has been found that the stationary state changes discontinuously as a function 
of the strength and the frequency of the external field. This is called the optical bistability [l9l - [2lj . Recently, a phase 
transition has been also predicted when the strength of interaction is driven periodically [22| . Furthermore, thanks 
to the recent realization of the atom-photon coupling in various kinds of experiments, many cooperative phenomena 
are studied extensively from experimental and theoretical views [2^. 

In the present paper, we study cooperative phenomena in a cavity due to the strong interaction and strong external 
field in a dissipative environment. As a model of a cavity, we adopt the Tavis-Cummings model 3] and the Dicke 
model j2^ . in which many two-level systems couple with a single quantized mode of the cavity photons. In order to 
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describe the system under a dissipative environment, we use a quantum master equation. In the hteratures [25l. [26}. 
the simple Lindblad equation is adopted. In those cases, the effect of the atom-photon couphng does not reflect on 
the dissipative term. Such a treatment is enough to study the system with relatively weak couplings. Indeed the 
optical bistability has been well studied with the simple Lindblad equation [2^, [2^. However, within this scheme, 
we cannot reproduce the ground state correctly as a stationary state in the case of no external field. In order to 
study overall structure of the phase diagram parameterized by the strength of the interaction and the intensity of the 
external field, we construct an extended master equation in which the effect of the interaction is taken into account. 
With this master equation, we reproduce the ground state correctly as a stationary state of the master equation. 

We obtain a stationary state by using the master equation and classify the stationary states into three phases 
by taking into account the symmetry of the time-evolution operators. The three phases are analyzed by the Fourier 
spectrum of the long-time behavior. We find the following three phases: "regularly oscillating phase" when the system 
is oscillating with the external field, "ordered phase" when the system has a spontaneous order, and "irregularly 
oscillating phase" when the system shows an irregular oscillation. 

We obtain two results: first, we obtain the phase diagram of the Tavis-Cummings model and the Dicke model under 
the AC external field. We find that the two models show a qualitatively different behavior when the interaction and 
the AC external field are strong and that the rotating wave approximation is not appropriate in this region. Second, 
we find a novel symmetry breaking phenomenon which is induced by the AC external field in the case of the Dicke 
model. 

The paper consists of the following sections. In Sec. II, the model of the system is explained. In Sec. Ill, a 
master equation which takes into account the atom-photon coupling is derived. In Sec. IV, we derive the equations 
of motions with mean-field (MF) treatment, and we show that ground states are realized in our approach in the case 
of no external field. In Sec. V, we study the phase diagram in the case of the Tavis-Cummings model. In Sec. VI, we 
study the phase diagram in the case of the Dicke model. In Sec. VII, the summary and discussion are given. 



II. SYSTEM HAMILTONIAN 



We study the system of photons and atoms in a cavity with an external field given by 

n (t) ^ ^ 

= cjpfl^a + wa^f + 5 X! { («^^ + °^^^) + ^ ("-^^t + aS-) } + 2^ {a) + a) cos {uj,t) . (1) 
1=1 1=1 

The first term represents the single quantized mode of the cavity photons with the resonance frequency uj-p which 
depends on the form of the cavity. The annihilation and creation operators for the cavity photon are denoted by a 
and a\ respectively. These operators satisfy the bosonic commutation relation [a, a^] = 1. In the second term, we 
describe many atoms as an ensemble of two-level systems with the energy gap Wa which are described by spin one-half 
operators. The spin operator Si(i — 1, • • • ,N) represents the ith two-level system with the commutation relation: 
[St,Sl] = 2Si. The third term gives the interaction between photons and atoms, and g is its coupling constant 
between them. The last term describes the AC external field with the frequency Wc and the amplitude ^. 

We study the time evolutions of the following quantities: photon field in the cavity (a), the number of excited atoms 
which is described by the z-component of the total magnetization (X]i=i ^i)^ and the dipole moments of atoms which 
are described by the transverse-component of spin operator {J^iLi ^t)- present study, we study ^-dependence 

for the case of Wa = = Wo, though lUc dependence on photon emission is an interesting problem and has been 
studied extensively in experiment [2l| . 

The Tavis-Cummings model is given by applying the rotating wave approximation to the Dicke model, which 
corresponds to the case x = while the Dicke model is defined by the case x = 1. In the two models, the nature of 
symmetry is different, the Tavis-Cummings model has U (!) symmetry where the ordered components (a) = |(a)|e''^ 
and J2i=i i^t) — ~ISi=i {^t) /N\e~^'^ can have any phase (j). On the other hand. The Dicke model has Z2 
symmetry where the ordered components are given by (a + a)) and Ylif=i {^t + ^i) 1^- This difference is important 
when we study the symmetry breaking phenomena under the AC external field as we will see in Sees. IVTl and IVIH 

In the present work, we study the macroscopic properties in the thermodynamic limit: N ^ 00 and V — >■ 00 keeping 
N/V to be a constant. Here, F is a volume of the cavity. In this case, we expect that the expectation values of a^a, 
SiLi '^i^ and X^ili '^t of 0{N). Consequently, the expectation value of the term '^f^iiaS^ + a'^ S~ ) + xia^ + 
aS^) is of 0{N'\/N). To obtain the meaningful result in the thermodynamic limit, we introduce scaled parameters. 

First, we introduce a scaled interaction between photons and atoms as 



(2) 
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The coupling constant g in Hamiltonian ([T]) is given by 



(down| e xj'-* • e |up) 



(3) 



where x^*'' denotes the position of the jth electron in the ith atom, and |down) and |up) express the eigenstates of 
Sf = —1/2 and = 1/2, respectively. We assume that the value of g does not depend on atoms. Here, e is the 
elementary charge and e is the polarization vector of the cavity photon. Because g is inversely proportional to the 
square root of the volume of the cavity V ^ the scaled value g is of 0(1) in the thermodynamic limit. With this scaling, 
the interaction term g{a + a^) Y^f^i{Sf + S~) scales as 0{N). 
Next, we also rescale the external field as 



N 



(4) 



Here, we assume that the amplitude of the external field increases proportional to v N. Then, the Hamiltonian is 
expressed with the new parameters as 



Hs{t) 



N 



: UJpG^a + J2 ^a^f + ^ H [{""^t + + X {a^S+ + aS^)] + 2VN( {a^ + a) cos {uj,t) . (5) 



N 



1=1 " " ■ i=l 

In the present paper, we study the nature of this Hamiltonian in the thermodynamic limit. 



III. MASTER EQUATION 



A. Stationary state in dissipative environment 



We investigate the time evolution of the system in a dissipative environment. The dissipation comes from the 
interactions between the system and the thermal baths. The master equation has the form 



(6) 



where r[ps] is a linear operator on pg for the relaxation process. 

After a long time evolution, the system reaches to a kind of stationary state at which the external field and the 
dissipation are balanced. In literatures [1^ [2^ , the following simple Lindblad equation 



N 



N 

T[ps] = -K {[a\aps] + [psa\a\) - 7,, ^ {[S+ , ps] + [psS+ , S']) - 7. ^ ([^f , Sfps] + [psS^, 5f]) , (7) 

i=i 



was adopted. In this case T[ps] represents the photon damping, the longitudinal damping of atoms 7// = {jxy + 27^), 
and the phase dephasing -f^ = 2'-fxy In this form, the effect of the atom-photon coupling in the system Hamiltonian 
does not reflect on the damping process. This effect was taken into account in some studies for the case of a single 
atom interacting with a cavity mode [27l[28j. while it has not been considered in a multiatom case. In the multiatom 
case, it is more important to take into account this effect because this effect changes the nature of phases as we will 
see in the present paper. 

In the present model ([51), it is known that when g is 0{y/uJ^[uJ^), the Dicke transition occurs in equilibrium = 0) 
by changing the coupling strength g [l3, [13 • When g is greater than a critical value 0{^ujaUJp), the photon field 
and dipole moments of atoms appear spontaneously. This ordered region is called "superradiant phase" . Because the 
relaxation terms in Eq. ^ always cause the reduction of photons and the relaxation of excited atoms, the master 
equation cannot correctly reproduce the ordered phase as a stationary state as we will see in Sec. IV. C. In the case 
of the Tavis-Cummings model, this simple Lindblad equation fails to have an ordered state at all. In the case of the 
Dicke model, it gives a phase transition, but fails to give the correct values of the order parameter. That is why it 
is necessary to obtain the master equation by projecting out the degrees of the freedom of the thermal baths, which 
incorporates effect of the atom-photon coupling into the relaxation terms. 
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B. General scheme for the master equation 



We have derived the master equation by tracing out the degrees of the thermal baths. For simphcity, we consider 
the thermal baths at zero temperature for photons and atoms, respectively. As for the coupling between atoms and 
a dissipative environment, we may adopt several types of interactions. In reality, because the range of the dissipative 
environment is spatially, a uniform coupling with all atoms takes place. We consider a thermal bath which models 
this effect, and call it "global-coupling bath" . When there only exists a global-coupling bath, the system conserves 

the total spin, — J2^=i{St '^^ + ^t)/'^ + {^tY ■ On the other hand, if we take into account the thermal bath 
coupling with each atom independently, the total spin is no longer conserved. In this case, we call the bath "local- 
coupling bath" . In order to study the relaxation phenomena, we need to take into account the local-coupling bath. 
In principle we have to contain both of the global and local coupling baths. In Appendix. [Cl we study dependence of 
the stationary states on the type of dissipative environment. If we include only the global-coupling bath, the system 
has an additional conserved quantity and the qualitative nature of the stationary states changes. In contrast, we find 
that the qualitative nature of the stationary states is similar as long as a contribution of the local-coupling bath is 
included. Therefore, in the present paper we adopt only the local-coupling bath as in previous works |25l. |26|. which 
represents the latter cases. 

Each thermal bath consists of free bosons and the Hamiltonian of thermal baths is given by 

N 

where An (wq) is the boson operator of each bath. We use the subscript 71 = for the photon and n = 1, . . . , TV for 
the nth atom. The interaction Hamiltonian is given by 

N 

■Hl^nXY, XnYn, (9) 
n=0 

where 

Xo = a + a\ X^^Sr + S+ (z = 1, 2, . . . , TV), (10) 



r„ = ^ kn {uJa) AI^ (uJa) + K i^a) ^ i'^a) , 
w„(>0) 



(11) 



and A is a coefficient to describe the strength of the coupling. In the present study, for simplicity, we consider the 
situation where the relaxation is only due to a radiative process, and thus we do not take into account the dissipation 
through . The quantum master equation of the system is derived in the standard way [29l - [3l| . As a result, we 
obtain 

dps jt) i 



where 



du 



kn (w) 



X (e-'"" [Xn, Xn (t, t-u)ps {t)] - e'"" [X„, ps {t) X„ (t, t - u)]) 



(12) 



X„(i,r) =Texp(^y dr' -Cs {r') ] X^, 
2 2 

fc„(w) = 2^ \kn{uJa)\ Su:^,i, 



(13) 
(14) 



Here, we define the Liouville operator as iCs (t) ■ = —i[Hs (t) , •]. 

In the derivation, we consider the Van Hove limit (weak-coupling limit) [32*1 where A — ^ keeping tX"^ to be constant. 
In this limit, the relaxation time of the system into the stationary state O (l/A^) is much larger than the typical 
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correlation time between the system and the thermal bath, and it is justified to neglect the contribution from the 
higher terms of Hi than the second order and to apply the Markov approximation j33!] . We also i gno re the term which 
represents the influence of initial correlations between the system and baths, which is shown in [3J| in the Van Hove 
limit to give the same result as one which a general initial state gives. 

In addition, in order to approximate the master equation into tractable form, we consider the case when correlation 
time is much smaller than the period of the external field. In this case, we may write 



Xn [tjt — m) ~ exp ( —-Hs (t) u I Xn exp ( -V-s [t) u 



(15) 

Here we regard the external field as a constant during the time evolution from t ~ u to t. With this approximation, 
the master equation is given by 

-^^-^[-Hs it), PS it)] 



A 



N 
n=0 



du 



duj 



Xn, exp --"Hs (t) u X„ exp -Tis (<) u ps (t) 



Xn,ps (t) exp -t'^s (t) u Xn exp -Hs [t) u 



iCt (t) PS- 



(16) 



It is noted that the master equation is applicable regardless of the number of atoms and the scaling of the coupling 
constant. If we ignore effects of the interaction and the external field in the Hamiltonian, we arrive at the form ([7]) 
(Appendix El). 



C. The stationary state and the averaged values of physical quantities 



In the present paper we study dynamics evolved by Eq. (|16p. Because the system is driven by the external field, the 
meaning of the stationary state is obscure. However, when we consider the symmetry of the time evolution operator 
CT{t), it is possible to define the stationary states and obtain the relation on the averaged values over the period of 
the external field. 

For the definition of the stationary state, we may use the Floquet eigenstates because the time evolution operator 
iCT (t) is periodic in time, 

iCr {t + To) = iCr (t) (t, = |^ j . (17) 
The eigenstate of Floquet operator J' is given by 

T(j)F (0) EE Texp ^ " dt'Cr (i')j '/'f (0) = e-"^'^'=(/)F (0) , (18) 

where 0f (0) and ep are the Floquet eigenfunction and quasi-eigenenergy, respectively. We expand ps (0) with the 
Floquet eigenfunctions. 

Ps(O) = 5]cf0f(O), (19) 

F 

where cp is a coefficient. We define the time evolution of the Floquet eigenfunction (/)f (i) as 

(t) = Texp (^i Ct {t')dt'^ 0F (0) . (20) 

In this definition, for all t, 

PS (t) = ^F'^F (^) ' (21) 
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and 



(/-F (t + Tc) -e-''^^^"</.F (t) 



(22) 



If the imaginary part of ep is zero, the corresponding state (pF (t) remains in the long-time hmit. We therefore define 
the stationary state ps,ss (t), 



PS,ss (t) = ^F-^F ' S = {F I Imep = 0} 



(23) 



Fes 



where S is a set of eigenfunctions with Im ep = 0. 

To discuss the relation on averaged physical values over the period of the external field, we introduce the unitary 
operator 



U = exp 



N \1 

/ . 



(24) 



which changes the sign of the operators a, , , and such that W all = —a. Under this unitary transformation, 
we have the relation 



(25) 



because of the following facts. In the Hamiltonian ([5]), all the terms except the last one are invariant with the 
transformation of U ^ while a + of the last term is changed to — (a + a^). Under time translation t ^ t + To/2, 
cos {ujct) changes the sign and the last term is also invariant. Thus, Eq. holds. With this property, we can show 
that U(j)F{Tc/2)W is an eigenfunction of T with the same quasi-eigenenergy of 0f(O): Because of Eq. ([TS]). 



3Tr 



(26) 



Now we consider averaged values (/ (a, , , S-))t in the stationary states. Here (. . . )t means Tr . . . ps,ss (t), and 
/ is an arbitrary function of operators. 

If the stationary state is unique, the set S consists of only one element Fq and the stationary state is given by 



Ps.ssit) = 0Fo(i)- 



(27) 



In this case, efq must be zero because Trpg^ss (t) — 1 must be satisfied. We find from Eq. (|26p and the uniqueness of 
the stationary state that 



0Fo(O) = C/0Fo 

In this case, we conclude 

(/(a, at, 5±, 5f ))i = Tr/(a, a\S^, 5f )ps,ss(0 

= Tr/(a,a^5'f ,S'f)Texp (^i dt'Crit' 

= Tr/(a,a^ S^,S^)Texp dt'Crit' 

dt'CTit')UFo 

= Tr/(-a, -a\ -^±,5^^3,38 [t + 
= {fi-a,-a\-S^,S-))t+Tj2- 



(28) 



MO) 



(29) 
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Averaging over the period of the external field, we have the relation 

(/ (a,at,5±,5f)), = (/(-a,-at,_5±,5,f))t, (30) 
where = • • • dt' jT^- For example, in the case that / (a, a} , Sf^ Sf) = a, 

Wt=V^t = -Wt, (31) 

and thus 

Wt = 0. (32) 

In a similar way, we can show 

WYt = 0. (33) 

Thus we find that averaged values of odd moments are zero if the stationary state is unique. 

Depending on the structure of the set S, there are several types of stationary states. In the case that the set S consists 
of several elements, the stationary state is given by a linear combination of those modes and it looks irregular. In the 
thermodynamic limit, the system may exhibit a phase transition. In this case, the set of S is degenerate reflecting 
the structure of the order parameter, and one of them is chosen as the symmetry breaking. For example, in the 
superradiant phase the spontaneous value of photons (a) and atom X^iLi i^t) appear, and the Floquet state in the 
set S corresponding to the symmetry broken state is realized. For the optical bistability which is shown in Sec. IVI[ 
two Floquet states exist in the set S in the thermodynamic limit. In Sec.|Vl we classify the stationary states according 
to this picture. 



IV. MEAN-FIELD TREATMENT 



In order to solve the time evolution, we need all the eigenvalues and eigenstates of the system [33]. It can be done 
only for systems with small degrees of freedom [13, [M HI, 113 ; It is difficult in the present case because the model of 



the cavity consists of many degrees of the freedom. As a solution of the difficulty, we focus on the nature of a cavity 
that all the atoms coherently interact with a single mode of photons, which causes an effective long-range interaction 
between atoms and study the time evolution of the system by a mean- field (MF) treatment. The MF treatment is 
rigorous in this case in the thermodynamic limit [s^. In the MF treatment, we assume that, for all the t, the density 
matrix is given by a product of density matrices of photons and atoms: 



PS 



it)^Ppit)®P^^Ht)^---®P^fHt), (34) 



where Pp{t) is the density matrix for photons and pi\t)[i = 1, • • • , N) is a 2 x 2 matrix for the density matrix of the 
ith atom. Here, we assume that all the atoms are in the same state Pait). In this case the time evolutions of pp (t) 
and Pa (t) are described by the Hartree-type equations. As a result, we can reduce the many-body problem to the 
one-body problem with MF parameters. 



A. The master equation for photons 



First we derive the master equation for photons. The equation is given by projecting out the degrees of freedom of 
atoms. We define the Hamiltonian for photons "Hp (t) by replacing the operators of atoms into the average values. 



^ 2—1 

+ 2Vn^ (a^ + a) cos (wct) + C {t) . 
By introducing a dressed bosonic operator 

N 

-4^ E ((^D + xiSt)) + 2v^ecos (a.ci) 



6 = a H 



(35) 



(36) 
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which satisfies the commutation relation [6, 6^] =1, this Hamiltonian can be rewritten as 

Hp (t) = hujpb^b + C' (i), 



(37) 



where C {t) and C" (i) are not the operators, and they do not contribute to the resuh. Thus, we redefine the 
Hamihonian as 



Tiu = hoj^jb^b. 



(38) 



This shows that the system is regarded to consist of the dressed photon field which incorporates the mean field of 
atoms. We derive the master equation of the cavity photon using Eq. (|16p . 



dt 



du I duj 



fco (w) 



{e-^"" [Xo, exp (-zHp(i)w) exp (zHp(t)M) Pp{t)] 
- e'"" [Xo, Pp{t) exp {-inp{t)u) Xo exp (iHp(t)M)] }, 
:-^[Hp(^),Pp(^)] 



duj 



ko{uj) 



X n{6{Lu-ujp) ([a + a^6pp(^)] - [a + , ppitp]) 
+ S{uj + ojp) ([a + a', b^Pp{t)] - [a + a\pp{t)b])} 

- i{p (^—) [a + a\bpp{t) + pp{tp] + P [a + a\b^ pp{t) + pp{t)b\ 

V — ujp J \uj + ujp J 

+ P 1=(1 + x) Y.^{S7) + {St)) + 4V]vecos(c.ei) [a + a\ Pp{t)\ }. 

Here, |fco(cLi)p = |^o('^a)pi5cjQ.Lj- Neglecting the contribution of the Cauchy principal value, we have 

^PM^_,^Up{t),pp{t)]~K{[a + a\bpp{t)\ - [a + at,pp(t)6t]), (40) 



(39) 



where we introduce a parameter k = ttA^ 



B. The master equation for each atom 

Next, we derive the master equation for each atom. Like the case of photon, the equation is given by projecting 
out the degrees of freedom of photons and other atoms. The Hamiltonian for each atom 'Kz. is given by 

Ha = huj^Sl + ^ [((a) + St + ((a^) + ^'l • (41) 



This Hamiltonian is diagonalized by a unitary operator [/a, 



where 



' + |^|(a)+x(at)|^ (43) 
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FIG. 1: (Color online) (a)Dependence of the ordered component |q| = |(a) /\/iv| in the stationary states on the interaction 
strength g in the Tavis-Cummings model with Wp = a;a = 1. Bullets show results of the present treatment with k = 7 = 0.1 
which agree with the exact results given by the red curve: a = ^ \/ " ^1/ 9^- Diamonds show the results of the treatment 
with simple Lindblad equation ((7]) with n = ^ = 0.1, where no ordered component appears at all. (b)Dependence of the ordered 
component |a| in stationary states on the interaction strength g in the Dicke model with ujp = uj^, = 1. Bullets show results of 
the present treatment with k = 7 = 0.1 which agree with the exact results given by the red curve: a — Ag'^ /uj'^ — i^l/ (4g^). 
Diamonds, squares, and triangles give the data obtained by the treatment with simple Lindblad equation ((Tj) with k = 7 = 
0.1, 0.01, and 0.001, respectively. They converge to a limiting value as k and 7 decrease, but the value does not agree with the 
exact value. 



and J7a is given by 



1 



^^ + ^l(«>+x(«t)|^ \^ 



m+x{a)) 
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with 



(44) 



(45) 



Here, we define |+) as eigenstate of a and |— ) as that of —a. As was seen in the case of photon field, the system is 
regarded as consisting of a dressed atom, Si{i — 1, ■ ■ ■ , N), 



(46) 



which incorporates the mean field of photon field. The master equation of the atom is given by 



dt 



i[Ha(i),Pa(t)] 

7{ (-1 Ul {S- + S+) C/a 1+) [S+ + Sr,S~.p,it) 



{+\UliSr + S+)U.\-) 



(47) 



Here, we introduced 7 = ttA^ \ki {2a)\^ as a constant because the energy scale of the bath is much larger than that of 
the interaction between the system and the bath. Like the case of the cavity photon, the relaxation term is related 
to the dressed atom Si. 



C. New aspects of the present master equation 



We derived the new types of master equations for photons Eq. PO]) and each atom Eq. (|47p . Here we compare the 
present treatment and the treatment with the simple Lindblad equation ^ with 7^ = 0. In the treatment with the 
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simple Lindblad equation, the effect of dissipation for (a) is described by 



= —Kia) 



(48) 



This means that the relaxation term causes simple reduction of the photon field at the constant rate k. 
the effect of dissipation for (Sf) is given by 



J^iSD diss 



= -jTvS^ {[S+,Srps] + [psS+,Sr]) 
= -2j ({Sn + ^ 



(49) 



The number of the excited atoms relaxes to zero, which is expressed as (S'f) — —1/2, at the constant rate 27. These 
damping terms are verified when the cavity photons and atoms are independent. But in the composite system of 
photons and atoms with large coupling constant g, they are not applicable because the system dissipates into the 
different state from (a) = and (Sf) = —1/2. 

In our treatment, the dissipation terms make the dressed photon field b and dressed atom Si relax. It is the most 
important point that the relaxation is related to these dressed values, which relaxes the system into the genuine 
ground state. The effect of dissipation for the dressed photon field b is given by 



-KTTb{[a + a^,bpsit)] - [a + a\psit)b^) 
-K.((6)-(6t)). 



(50) 



This shows the relaxation of the dressed photon field. The effect of dissipation for the dressed atom S'f is given by 



d ,~ 
dt 



(^f )diss = -7Tr^f ((-I Ul {Sr + S+) |+) [s+ + Sr , ps{t)] + h.c.) 



(51) 



which describes the relaxation of an excited dressed atom. Here, the coupling constant 7 is also modulated by the 
mean field as 



7[s,(a>]^l {-\Ul{S- +S+)U,\+) p7- 



(52) 



With these dissipation terms, the system indeed relaxes to the true ground state of the Hamiltonian in the case of 

in the stationary states 



no external field (^ = 0). We show the dependences of the ordered component \a\ ~ (a) /vN 

on the interaction strength g of the Tavis-Cummings model and the Dicke model in Fig.[lja) and (b), respectively. As 
we see, the results of the present treatment (bullets) reproduce the equilibrium states (red curves). In the figures, we 
also give the results obtained by the treatment with the Lindblad equation ^ with various strength of the couplings, 
K = 7 = 0.1 (diamonds), 0.01 (squares), and 0.001 (triangles). In the case of the Tavis-Cummings model (Fig. [ija)), 
the equation (O does not show any ordered state. In the case of the Dicke model (Fig.[TIb)), the simplified equation ([7]) 
gives a phase transition, but fails to give the correct values of the order parameter. 

Thus, the present formalism of dissipation terms satisfies the minimal condition for the study of the strong inter- 
action cases. The terms are regarded to be the extension of the master equation with dressed modes. Therefore we 
use this formalism in the present study. 



D. The equations of motions 



We express the master equation with the scaled quantities, 
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The equations of motions for a, m"*", and are given by multiplying operator a, , and Sf respectively on the 
corresponding master equations Eq. (|40p and Eq. (|47p . respectively. As a result, we obtain a set of classical equations 
of motions: for the photon field, 



da 
'dt 



for the dipole moment of atoms 



i (wptt + g (to + x™^) + 2^ cos {ujct)) 
H (to + xm 



— (to + 



(54) 



2i I — ^ ('^* ^ TO^I 



|r72+g2|Q,_)_^Q,*|2| 



+ {^n^ ~ .9^ (" + XQ?*)^! |w (a* + X"^) + v'^m'^ + 9^ {a* + xo^f m | 



(55) 



and for the number of the excited atoms 



dt 



ig {{a + XQ^*) ^^t-^ ^ {<^* + X^^) } 
7 



+ g2 |q, ^ T^Q,* |2| 

{r;2 - g^{a* + xaf} i^w (a + xa*) (m+ + to") - g^ {a + xa* f (^m' " ^) + ("^' + I 
- (a + xa*)'} | V9 (a* + X") (™+ + "i") + t?' (^to^ + ^ " («* + X")' (^™' - ^) 



(56) 



In the region, g ^ Wa, g ^ Wp, and ^ ^ Wp, with the rotating wave approximation for the interaction term between 
the system and the bath, these equations are reduced into those obtained by the simple Lindblad equation ([T])- These 
equations seem to be complex, but the most important thing here is that the effect of the atom-photon coupling is 
incorporated into the relaxation terms. The nonlinearity of the equations increases with ga, which causes the Dicke 
transition and complicated phenomena in large ^. In the next section, we discuss the behavior of the stationary states. 

V. TIME-EVOLUTION AND TYPES OF STATIONARY STATES 

The stationary state in thermodynamic limit can be obtained by evolving the state by the equations of motion (l54l) . 
(f55| . and (|56|) for a long time instead of obtaining the floquet state ps.SS {t) explicitly. We concentrate on the case 
uja = ujp = u)c = uj which is set to be one and k = 7 which is set to be 0.1. We solve the equations by means of 
Runge-Kutta method. In the present study, we set time interval At — 0.00027r which is fine enough for the present 
simulation. We confirm that this value of At is enough small by changing At. We performed simulations with 
increments Ag = 0.02 and — 0.001. To observe the stationary state, we study the behavior in the large t region, 
e.g., {t = 98007r — IOOOOtt). We confirmed that the state reached a stationary state by changing the range of t. 

In order to classify the stationary states, we make use of the Fourier transformation to detect this symmetry. 



1 



IOOOOtt 



SOOvr igsoOTT 



dt'a {t') e''"^"*' {k = O.OIto | 



TO 



e Z,-210 < TO < 210). 



(57) 



According to Eq. ([30|) . if the stationary state is unique. 



/ (a, a* , , m^) — f {—a, —a* , , m^), 



(58) 
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and only the odd numbers of k in Eq. (j57l) are allowed: 

a (t) = Qffc exp (ikujct) ■ (59) 

fc— (odd numbers) 

On the other hand, if the symmetry is broken, then the even numbers are allowed, too: 

a (t) = ttfe exp (ikujct) ■ (60) 

fc— (integers) 

We call the static component ao "spontaneous order of photon" . 

When the symmetry is not broken, ak has delta peaks at k = ±1. This phase represents that the system is only 
driven by external field and we call this "regularly oscillating phase" . If the symmetry is broken, we observe the delta 
peaks at fc = and k = ±2 in addition to those at fc = ±1. In this phase, the spontaneous order of photon ao appears 
and we call this phase "ordered phase". Besides these phases, there are regions which do not belong to the above 
two phases. In these regions, the Fourier spectrum has delta peaks at non-integer fc, or even shows broad peaks. We 
call this "irregularly oscillating phase" . For each Fourier component fc, we set the threshold for the existence of delta 
peak to be \ak\^ = 1.0 x lO""^. 

In experiments, we can measure the photon emission from the cavity which represents the energy transfer per unit 
time from the cavity to the thermal photon bath. It depends on the dressed value (see Appendix |B]) given by 

Ip = iKUjpilmbf. (61) 

It is possible to classify three phases by the photon emission. In the regularly oscillating phase, the Fourier coefficients 
have delta peaks at fc = and fc = ±2. In the ordered phase, they have delta peaks at fc = ±1 in addition. In the 
irregularly oscillating phase, they also have delta peaks at non-integer fc, or even shows broad peaks. 

VI. PHASE DIAGRAM OF THE STATIONARY STATES OF TAVIS-CUMMINGS MODEL 

In this section, we study the phase diagram in the Tavis-Cummings model driven by the AC external field in the 
coordinate of the strength of interaction g and the intensity of external field ^. Here, we show two results: First, we 
confirm the optical bistability to appear in the region g ~ a; and ^ w, which has been studied in the region where 
(jf <C w and ^ <C w (20l . [2ll [26j . Second, we observe the irregular oscillating phase in the region where g and ^ are 
strong. 

We study the stationary state when ^ is increased and when ^ is decreased at a fixed value of g, respectively. When 
we change ^, as the initial state we use the stationary state of a, m^, and to"*" for the previous value of ^. We divide 
the phase diagram into three phases: regularly oscillating phase, ordered phase, and irregularly oscillationg phase. 

From now, we will explain the phase diagram in Fig. [2] in the following three subsections. In the first section fVIAl 
we will study the region where the optical bistability appears. In this region, the system is in the regularly oscillating 
phase. In the second section fVIBl we will explain the strongly interacting and weakly excited region. In this region, 
the effect of interaction between atoms through photon field is dominant. Therefore, we find an ordered phase due to 
the strong interaction g (Dicke transition). Finally, in the third section IVI C[ we will show the behavior in the highly 
interacting and highly excited region. Here, we find the irregularly oscillating phase. 

A. Optical Bistability 

In the region where g and ^ are weak {g^ < 0.07), we find only the regularly oscillating phase. In this phase, the 
photon amplitude a oscillates around the origin and the dipole moment of atoms oscillates around the origin in the 
(m^, TO^) plane. 

In Fig. |3] (a) , we show the ^-dependence of the photon emission /p averaged over the period of the AC external field 
in a case of very weak interaction g = 0.12. In this figure, Ip smoothly changes with increasing and decreasing ^. 

In Fig. [3] (b), we also show the ^-dependence of Ip for g = 0.4, where we confirm the optical bistability to appear; 
We find a discontinuous transition in each of the ^-increasing process and in the ^-decreasing process, and the points 
where the transitions appear are different. Namely, we find a hysteresis phenomenon. As for m^, jumps up from 
—0.2 to around zero in the ^-increasing process and jumps back in the ^-decreasing process at the different points from 
that in the ^-increasing process. It is the transition from a regularly oscillating phase to another regularly oscillating 
phase; the static order parameters do not appear. The optical bistability is a kind of non-equilibrium phase transitions 
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Tavis-Cummings model ( x=0-0 ) 




FIG. 2: (Color online) The phase diagram of the Tavis-Cummings model is parameterized by the interaction strength g and 
the intensity of the external field ^. This phase diagram shows the steady states in the processes of increasing or decreasing ^ 
at the fixed value of g for the case of a;a = i^p = i^c = 1. (i) The black bullets and red bullets denote points of discontinuous 
changes of the state in the ^-increasing and ^-decreasing processes, repectively (see Sec. IVI A|) . The optical bistability takes 
place between them, (ii) The regularly oscillating phase is denoted by blank space (see Sec. IVI (iii) The ordered phases in 
the ^-increasing and ^-decreasing processes are denoted by black thin lines and a red dotted lines, respectively (see Sec. IVI B|) . 
(iv) The irregularly oscillating phases in the ^-increasing and ^-decreasing processes is denoted by black bold lines and red bold 
lines, respectively (see Sec. lVIC)l . 



(a) 



(b) 



g=o.i2 




g=o.4 







. ■* 

^ i 


.'•'/ 







0.2 0.4 0.6 5 



0.14 0.16 0.18 5 



FIG. 3: (Color online) Dependence of the energy dissipation /p on ^ for g — 0.12 (a) and for g — 0.4. 



originated in the relaxation process. The external field continuously pumps photons vifhich are to be relaxed through 
two processes, one of them is via a direct photon relaxation processes, and the other is via the atoms. As we find in 
Eqs. (j54p and (j55p . when g is small compared to oja and Wp, the relaxation of photon via the direct photon relaxation 
is proportional to na, while the relaxation of photon via atoms is of order g^am^ /^. Therefore, if the external field is 
weak, both the relaxation processes contribute because w ~l/2, because the system is still near the ground state. 
When the external field is strongly applied, it is mainly through the photon emission because the energy emission via 
atoms is suppressed and approaches to zero. Between these cases the mechanism changes discontinuously. 
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B. Ground state phase transition due to the interaction g 



g=1.04 




0.4 \ 



FIG. 4: (Color online) Dependence of the ordered component |a| on ^ for g — 1.04. 



= ^WaWp — 1, the photon field and dipole moments of atoms appear 
In Fig. |4l we show the ^-dependence of the order parameter |a| for 



When g is larger than the critical value 
spontaneously as the Dicke transition (TtI. Tl8|. 
g = 1.04. We find that the ordered component of photon disappears continuously. When ^ is increased, the ordered 

component of atoms m+ also disappears continuously, and approaches to —0.5. In both processes of increasing 

and decreasing ^, we do not find hysteresis. In this case, the AC external field disturbs the order, which is plausible 
because the external field excites the system from the ground state. 

In the phase diagram, we draw the black thin lines in the ^-increasing processes and the red dotted lines in the 
^-decreasing processes, where the ordered component appears. 



C. Highly excited and highly interacting region 

In the region where g and ^ are strong (g^ > 0.07), we find the irregularly oscillating phase. In this phase, wc do 
not observe the order parameter. In Fig. \E\ we show the time evolution of Ip for g) = (0.37, 0.86). 

In the phase diagram (Fig. [5]), we draw black bold lines in the ^-increasing process and red bold lines in the ^- 
decreasing process, where the irregularly oscillating phase appears. On the phase boundary between the irregularly 
oscillating phase and the regularly oscillating phase, g is almost proportional to ^. The detailed analysis of this 
property will be studied in the future work. 



g=0.86 5=0.37 




30500 



30600 



30700 30800 t 



FIG. 5: (Color online) Time evolution of Ip in the region of irregularly oscillating at ^ = 0.37 for g — 0.86. 
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VII. PHASE DIAGRAM OF THE STATIONARY STATES OF THE DICKE MODEL 



g 



Dicke model ( x=l-0 ) 



0.6 



0.4 



0.2 



0.2 



0.4 



0.6^ 



FIG. 6: (Color online) The phase diagram of the Dicke model is parameterized by the interaction strength g and the intensity 
of the external field ^. This phase diagram shows the phase of the stationary state in the process of increasing and decreasing ^ 
at each fixed value of g for the case of toa = tOp — toa = 1- (i)The black bullets and red bullets denote points where the stationary 
state changes discontinuously from the regularly oscillating phase to another regularly oscillating phase in the increasing and 
decreasing ^ processes, respectively. The red squares observed around {^,g) = (0.13,0.54) show the discontinuous jump from 
the ordered phase to another ordered phase in the ^-decreasing process. (ii)The regularly oscillating phase is denoted by the 
blank space. (iii)The ordered phase in the ^-increasing and ^-decreasing processes are denoted by black thin lines and red 
dotted lines, respectively (see Sec. IVII"X| . (iv)The irregularly oscillating phase in the increasing and decreasing ^ processes are 
denoted by black bold lines and red bold lines, respectively (see Sec. IVII"^ . 



Next, we will study the phase diagram in the Dicke model driven by the AC external field characterized by the 
strength of interaction g and the intensity of external field ^. In this section, we show the two following results: First, 
the rotating wave approximation is justified when g and ^ are weak, while it is not justified when g and ^ are strong. 
Second, instead of the irregularly oscillating phase in the Tavis-Cummings model, we find a novel symmetry breaking 
phenomenon due to the AC external field when g and ^ are strong. 

We also divide the phase diagram into three phases: regularly oscillating phase, ordered phase, and irregularly 
oscillating phase as we saw in the case of the Dicke model. In Fig. [HI we give the phase diagram in the same notation 
of Fig. H 

From now we will explain the phase diagram. When the atom-photon coupling and the intensity of the AC external 
field are small{g^ < 0.06), we observe a similar feature in the phase diagram to the case of the Tavis-Cummings model; 
The region belongs to the regularly oscillating phase. When the interaction is stronger, we find the optical bistability. 
It shows that the rotating wave approximation is justified and that the effect of irrotation term qualitatively does not 
appear in this region. 

We also find the similar behavior in the region where the interaction is strong and the external field is weak. When 

the interaction exceeds = the spontaneous order parameters \a\ and m'^ appear [l3,[3|- When ^ is 

increased, the order parameters decrease gradually and disappear continuously to zero as we saw in the case of the 
Tavis-Cummings model. 



A. Ordered state of driven component in the highly excited and highly interacting region 



In contrast to the Tavis-Cummings model, the ordered state of driven component appears in the highly excited 
and highly interacting region. In Fig. [7] (a) and (b), we show the ^-dependences of the order parameters |a| and 



16 



(a) 



g=0.34 



(b) 



m*| g=0.34 





0.2 0.4 



0.2 0.4 6^ 



FIG. 7: (Color online) Dependence of quantities on ^ in the region 0.1 < (, < 0.6 for g — 0.34. (a) The energy dissipation \a\ 
on ^, (b) the dipole moment \m+\ on ^. 



m , g=0.40, 5=0.41 
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FIG. 8: (Color online) Time dependences of (black bold line) and (red dotted line). The period of the system is twice 
as long as that of the AC external field. 



Order parameters appear and disappear repeatedly and |a| and m+ behave in the same way. In the 

increasing process, order parameters appear discontinuously and disappear continuously. In the ^-decreasing process, 
they exhibit the same way. In the both processes, the continuous transitions appear at the same value of ^ and the 
discontinuous transitions appear at the different value of ^. From the argument given in Sec. |Vl the finite order 
parameters indicate the spontaneous symmetry breaking. In the region photon field a and the dipole moment of 

atoms oscillate around the finite value. Here, it should be noted that the value of m+ exceeds the bound around 
^ = 0.5, which should be attributed to approximations in the derivation of the master equation. This problem will be 
studied elsewhere. 

Next, we show the behavior of the irregularly oscillating phase. As we saw in Sec. |Vl "irregular" means that the 
system does not have the same period of the external field. In Fig. [51 we show the time dependences on (black 
bold line) and to^ (red dotted line) for g — 0.4 and ^ — 0.41. We find that the period is twice as long as that of the 
AC external field, and is always negative, which shows the symmetry breaking. 

In the phase diagram, the repeated belt-like structure appears which indicates that the AC external field makes 
the order parameter align at the discrete values of ^. On the phase boundary, g is almost inversely proportional to ^, 
which will be studied in the future work. 



VIII. SUMMARY AND DISCUSSION 



In the present paper, we studied the cooperative phenomena of photons and atoms under the AC external field. 
In order to study the stationary state, we used the formalism of a master equation for open systems. The Lindblad 
equation ([7]) is not appropriate when the atom-photon coupling is strong compared to the energy of photons and 
atoms. We derived the master equation which incorporates the atom-photon coupling into the dissipation terms with 
the use of the mean-field treatment under an approximation that the correlation time of the thermal bath is regarded 



17 



to be short. The dissipation terms of the present master equation are related to dressed quantities; The dissipation 
terms relax the dressed photon b which incorporates the mean field produced by atoms and relax the dressed atom 
which incorporates the mean field produced by photons. We studied the stationary state without the external field 
= 0) and confirmed that the dissipation terms reproduce the Dicke transition correctly. 

We showed the following results. We presented the phase diagrams parameterized by the strength of the atom- 
photon coupling g and the intensity of the AC external field ^ for the Tavis-Cummings model and the Dicke model 
under the AC external field. For small values of the product g^, we found the similar behavior in the both phase 
diagrams. For small the system belongs to the regularly oscillating phase and we observed no singular behavior. 
When g increases a little, the optical bistability appears while changing the value of ^. When g exceeds the critical 
value, the spontaneous order parameters appear as the Dicke transition. In the region where the external field is weak, 
the AC external field disturbs the order. For large values of the product 17^, we found the different behavior in the 
phase diagrams for the Tavis-Cummings model and the Dicke model and therefore the rotating wave approximation is 
not appropriate in this region. In the case of the Tavis-Cummings model, the irregularly oscillating phase appeared. 
On the phase boundary between the irregularly oscillating phase and the regularly oscillating phase, g was almost 
proportional to In the case of the Dicke model, the repeated belt-like structure of the ordered phase appeared 
instead of the irregularly oscillating phase, which describes a non-equilibrium phase transition with a spontaneous 
symmetry breaking. On the phase boundary between the ordered phase and the regularly oscillating phase, g was 
almost inversely proportional to ^. 

The detailed analysis on the phase boundary which is approximately given by a line g/^— constant in the Tavis- 
Cummings model and by a line given by a line g^— constant in the Dicke model and the mechanism how the ordered 
states appear under the AC external field will be reported elsewhere. 
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Appendix A: The simple Lindblad equation 



In this section, we derive the simple Lindblad equation ([7]) from Eq. (|16l) . 



dps {t) i 



N 



n=0 



du 



duj 



kn (w) 



First, we drop the irrotational terms in the interaction Hamiltonian between the system and the thermal bath. It 
is expected not to be essential because the interaction strength between the bath and the system is weak. We then 
obtain 

dps jt) i ,^ 



kn (w) 



where 



Ag =a,Xf^ = a\ 
Xi —S^ , — s 



{i^l...N). 



(A2) 

(A3) 
(A4) 
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Next, we drop the atom-photon coupHng and the external field in the system Hamiltonian in the relaxation terms: 



N 



(A5) 



which is justified in the region g <C 5 <C Wp, and ^ <C Wp. If we neglect the contribution of the Cauchy principal 
value, we obtain the simple Lindblad equation (jT)) given by 



dps (t) 
dt 



N 



- [ns it) , PS (<)] - « ( [a^ aps (t)] + [ps (t) a\ a]) - jj^ i ^ ^^Ps (t)] + [ps (t) S+ , S']) , (A6) 



where k = ttA^ 



ko (wp) 



and 7 = ttA^ 



h (Wa) 



It is noted that we can not neglect the effect of the interaction term and the external field easily because Wa — We 
and iOp — ujc are important in the stationary driven state. In the resonant case where = i^a — '^p; the system indeed 
shows the optical bistability as a cooperative phenomena via atom-photon coupling. 



Appendix B: the evaluation of the intensity of the transmission light 

In this paper we classify the stationary states into three phases. It is possible to distinguish them by the fourier 
spectrum of the intensity of the transmission light from the cavity /p. In this Appendix, we derive the specific form 
of /p. We consider the time development of the energy of the system E{t), which is given by the master equation: 

- V/ du dc^A2|fc„H|2TrHsW(e"'""[^n,^„(t,i-")psW]"e'""[^«,Ps(i)^n(i,t-u)])(Bl) 

The first term represents the energy transfer by the external field and the second term represents the energy dissipation. 
The intensity of the transmission light Ip corresponds to the 71 = term of the energy dissipation, 

du / duX^\h{co)\^Ti-ns{t){e-'^''[Xo,Xo{t,t-u)ps{t)]-^^^[Xo,ps{t)Xo{t,t-u)]), (B2) 
Jo 

which represents the energy dissipation from the system to the photon bath. After making the mean-field treatment 
and regarding the external field as a constant during the time evolution from t — u to t, 



Nip = 



Nip = 



du dc^A2|fcoH|'Trphoton(e-*""Hp(0 [a + a^, e-'«p(*)" (a + a^) e^p^'Vp (t) 

-e'^''Hp{t)[a + at,pp(t)e-*^p(*)"(a + a^)e*^p(*)" 

KTrphoton'Hp (t) ( [a + a\bpp {t)] - [a + a\pp (t) ) , 
AKLOp (Im 



(B3) 



where k = 7rA^|fco (i^p) P and b ~ a+ gJ2iLi {{^i ) + x{S^)) /VN + 2-\/]Vf cos {uJct) /wp. In the last transformation, 
we approximate that photon state is coherent. With the scaled quantities, Ip is rewritten as 



Ip — 4:KUjp \ Im 



a + (to -|- 



(B4) 



Appendix C: The effect of the thermal bath coupling with all the atoms uniformly 



In this Appendix, we study the dependences of stationary states on the type of dissipative environment which 
consists of a global-coupling bath and local-coupling baths. The global-coupling bath is given by the following 
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FIG. 9: (Color online) Dependences of the order parameter of photons |q| on ^ in the ^-increasing process for the case of 
the Dicke model (x = !)■ We use black circles for 7 = 0.1 and 7t = (Fig. [Tjl, blue squares for 7 = 0.1 and 7t = 0.1, 
and red triangles for 7 = 0.1 and 7t = 0.2. We find the qualitatively similar behavior in these different types of dissipative 
environments. We also plot the points for 7 = and 7t = 0.1 by crosses and we find the different behavior from other cases. 



Hamiltonian 

N 

= ^a^TK)^TK) + A^(5+ + 5r) ^ (fc;(w„)A(c^,) + fcT(c<Jc.)A^K)), (CI) 

"q(>0) i=l cJq(>0) 

where the first term represents the free boson bath and the second term represents the interaction between all atoms 
in the system and the thermal bath. Here, A}^(uja) and At(wq.) are creation and annihilation operators of the thermal 
bath, respectively. In the same manner as in Sec IIVBI we obtain the contribution of the global-coupling bath to the 
equation motions of atoms by tracing out the degrees of the freedom of the thermal bath: for the dipole moment of 
atoms, 

d_ + ^ ^7T 2 

— [rf — g'^{a + XQ!*)^)(77^m^ 
and for the number of the excited atoms. 



2775(0; + xa*)™^) 

2r^g{a* +xa)m-)i (C3) 

where 7t = E<.^„(>o) \kT{uJa)\'^S^^^2a- 

In Fig. |9l we show the stationary states with various kinds of the coupling strength between the system and the 
thermal baths. Here, 7 and 7t denote the coupling strengths of the local-coupling bath and that of the global-coupling 
bath, respectively. We show the stationary states for 7 — 0.1 and 7t = by black circles, 7 = 0.1 and 7t = 0.1 by blue 
squares, 7 = 0.1 and 7t = 0.2 by red triangles, and 7 = and 7t — 0.1 by crosses. When a dissipative environment 
only consists of the global-coupling bath (crosses), the total spin is conserved, and thus we find a qualitative different 
behavior from the cases with non-zero 7. On the other hand, when the total spin is not conserved (7^0), we find a 
similar dependence. We also checked the effect of the global-coupling bath (7t) does not cause qualitative change in 
other parameter regions. 



2r]g{a + xa*)^^) 
— g'^ioi + xQ*)^m^ — 2rjg(a* + xq;)to^)] 



(C2) 



Q^^A^ss- (^2+^^2|^ + ^^*|2)2l^ i 

X [(jf' — g'^{a* + xQ^)^)('?^'7i^ — g'^{oe + xQ!*)^m^ — 
^ iv'^ ^ 5^(0^ + XQ;*)^)(^^"^^ ~ 5^(Q! + XQ!*)^m^ — 
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